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The Super Proton Synchrotron (SPS) covers one of the most interest- 
ing regions of the QCD phase diagram (T — /^b). On the one hand there 
are indications that the energy threshold for dcconfincmcnt is reached al- 
ready at low SPS energies. On the other hand theoretical calculations 
locate the QCD critical end-point at energies accessible at the SPS. In this 
paper the NA49 signatures of the onset of deconfinement are presented. 
Results are shown on pion production, the kaon to pion ratio, slopes of 
transverse mass spectra ("temperature"), and event-by-event particle ratio 
fluctuations versus collision energy, for central Pb + Pb interactions. Next 
we show possible indications of the critical point in event- by- event mean 
transverse momentum and multiplicity fluctuations. Finally we discuss the 
future ion program of the NA61 /SHINE experiment (energy scan with light 
ions). 
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1. Introduction 

It is a well established fact that matter exists in different states. For 
strongly interacting matter at least three are expected: normal nuclear mat- 
ter (liquid), hadron gas (HG), and a system of deconfined quarks and gluons 
(eventually the quark-gluon plasma). One of the most important goals of 
high-energy heavy-ion collisions is to establish the phase diagram of strongly 
interacting matter by finding the possible phase boundaries and critical 
points. In principle, we want to produce the quark-gluon plasma (QGP) 
and analyze its properties and the transition between QGP and HG. 

The phase diagram of strongly interacting matter (Fig. []JJ is most often 
presented in terms of temperature T and baryochemical potential ^b- It is 
believed that for large values of [Ib the phase transition is of the first order 
(gray band) and for low \xb values a rapid but continuous transition (cross- 
over). A critical point of second order (CP) separates those two regions. The 
open points in Fig. [1] are hypothetical locations reached by the high-density 
matter droplet after dissipation of the energy of the incident nucleons from 
where the evolution of the expanding and cooling fireball starts. The closed 
symbols are chemical freezout points [1] (chemical composition fixed, the end 
of inelastic interactions). The fireball then expands further until thermal 
(kinetic) freeze-out takes place (particle momenta are fixed). 
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The Super Proton Syn- 
chrotron (SPS) covers one of 
the most interesting regions of 
the QCD phase diagram (T — 
(1b)- On the one hand it 
is expected that the energy 
threshold for deconfinement is 
reached already at low SPS 
energies (see the open point 
hitting the transition line in 
Fig. [1]). On the other hand 
lattice QCD calculations locate 
the critical point of strongly in- 
teracting matter in the SPS en- 
ergy range [2]. 

2. NA49 experiment 

NA49 [3J, operating in 1994-2002, was one of the fixed target heavy- 
ion experiments at the CERN SPS. The main components of the detector 
were four large volume time projection chambers (TPC). The Vertex TPCs 
(VTPC-1 and VTPC-2), were located in the magnetic field of two super- 
conducting dipole magnets. Two other TPCs (MTPC-L and MTPC-R) 
were positioned downstream of the magnets symmetrically to the beam 
line. The NA49 TPCs allowed precise measurements of particle momenta 
p with a resolution of o-(p)/p 2 = (0.3 — 7) • 10~ 4 (GeV/c) -1 . A precise 
measurement of specific energy loss (dE/dx) in the region of relativistic rise 
was possible in the TPCs. Time of Flight (TOF) walls supplemented par- 
ticle identification close to mid-rapidity. The typical dE/dx resolution was 
a[dE j dx) I \dE / dx) « 0.04, the typical ToF resolution aiToF) w 60 ps, and 
invariant mass resolution (for identification of particles via decay topology) 
was <j{rrii nv ) ~ 5 MeV. The centrality of nuclear collisions was selected via 
measurement of the energy of the projectile spectator nucleons in the For- 
ward Calorimeter. The NA49 acceptance covers the forward hemisphere, 
but because of the symmetry of nucleus+nucleus (A + A) collisions this 
nevertheless allows to extract Air integrated multiplicities. 

3. NA49 evidence for the Onset of Deconfinement 

3.1. Kink, horn, step and SMES 




500 ^ 1000 



H b (MeV) 

Fig. 1. QCD phase diagram. 



In 2002 the NA49 experiment completed the SPS energy scan of central 
Pb + Pb collisions. This program was originally motivated by predictions 
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of the Statistical Model of the Early Stage (SMES) [I] assuming that the 
energy threshold for deconfinement (the lowest energy sufficient to create 
a partonic system) is located at low SPS energies. Several structures were 
expected within SMES: the kink in pion production (due to increased en- 
tropy production), the horn in the strangeness to entropy ratio, and the step 
in the inverse slope parameter of transverse momentum spectra (constant 
temperature and pressure in a mixed phase). 

Figure [2] shows production 
of charged pions (the total en- 
tropy is carried mainly by pi- 
onsQ) (vr) = 1.5((vr+) + (*-)) 
normalized to the number of 
wounded nucleons versus Fermi 
variable F (F « (stvjv) 1/4 )- 
In SMES, this ratio is propor- 
tional to the effective number 
of degrees of freedom (NDF) 
to the power of 1/4. For cen- 
tral A + A collisions (Pb + Pb 
for SPS or Au+Au for AGS and 
RHIC) a change of slope around 
30A GeV is visible (slope in A+ Fi S- 2 - Energy (F) dependence of the mean 
A increases from « 1 (AGS) to P ion multiplicity per wounded nucleon [5] in 
« 1.3 (top SPS+RHIC) - con- full phase space (Air). 
sistent with increase by a factor 

of 3 in NDF). Such an increase is not observed for p + p(p) reactions. The 
increase in NDF, when going from hadron gas to QGP, may be interpreted 
as a consequence of the activation of partonic degrees of freedom. 

Figure [3] (left) presents the (K + ) / (ir + ) ratio in full phase space (4-7r) 
versus energy. In SMES, the ratio is proportional to strangeness/entropy 
densities. Results for A + A are very different from the results for p + p 
and show a sharp peak in (K + ) / (ir + ) at 30A GeV. This peak is even more 
pronounced at mid-rapidity (see [5]). Recently, these intriguing NA49 re- 
sults on pion and kaon yields (at mid-rapidity) were confirmed at -s/snn = 
9.2 GeV and 19.6 GeV by the STAR experiment [BJ. The string hadronic 
models (curves in Fig. [3] (left)) do not reproduce the data. The hadron gas 
model (HGM) [7] using a parameterization of the energy dependence T c ^ 
and \ib based on fits to hadron yields produces a broad maximum in the 
(K + ) / (ir + ) ratio as a consequence of saturating T c h and decreasing /ig with 

1 In SMES the total entropy and the total strangeness are the same before and after 
hadronization (the entropy cannot decrease during the transition from QGP to hadron 
gas), therefore pions measure the early stage entropy. 
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Fig. 3. Left: Energy dependence of the (K + ) / (n + ) 
Right: E s measure of strangeness to entropy ratio versus energy in Air [5]. 



ratio in full phase space (47r). 



increasing energy (limiting temperature reached somewhere at the SPS). 
However, this version of the model overestimates the relative kaon yields 
from 30A GeV on. The latest extension of the HGM [8] yields and im- 
proves fit of the energy dependence (K + ) / (tt + ) by inclusion of controversial 
[9] a state and unmeasured resonances above 3 GeV in the model hadron 
spectrum (see conference slides or [8]). 

The measure which much better reflects the total strangeness to entropy 
ratio in the SPS energy range is E s = ({K} + (A))/(ir), proposed in Ref. [4|, 
and calculated from ir, K, and A yields in 4ir acceptance. The E s ratio 
can be directly and quantitatively compared to SMES predictions. Figure 
(right) shows a distinct peak in E s at 30A GeV. This behavior is described 
(predicted) only by the model assuming a phase transition (i.e. SMES), 
where the maximum, called 'horn', is the result of the decrease of strangeness 
carrier masses in the QGP (m s < m,A,K ...) an d the change in the number of 
degrees of freedom when reaching the deconfined state. 

Figure H] presents inverse slope parameters (T) of transverse mass spec- 
tra of positively and negatively charged kaons. For A + A data one can 
see a strong rise at AGS, plateau at SPS, and rise towards RHIC energies. 
Such structure is not observed for p + p collisions. The plateau is consistent 
with constant temperature and pressure in the mixed phase (latent heat) 
|10| . In fact, this structure strongly resembles the behavior of water, where 
a plateau can be observed in the temperature when heat is added. Models 
without phase transition do not reproduce the A + A data, but a hydro- 
dynamical model incorporating a deconfinement phase transition at SPS 



2 Transverse mass spectra were parametrized by dnj (mTdmr) = C ■ exp(—iriT /T); fits 
were done close to mid-rapidity. 



energies [TT] describes the results in Fig. 0] quite well. The step-like feature 
is also present in the energy dependence of mx — m of protons and pions 
(see conference slides or [5]). 
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Fig. 4. Energy dependence of the inverse slope parameter T of the transverse mass 
spectra of K + and K~ mesons [5]- 



3.2. Event-by-event particle ratio fluctuations 

Hadron ratios characterize the chemical composition of the fireball and 
are not affected by hadronic re-interaction when looking at conserved quan- 
tities such as baryon number, strangeness. We expect a change of particle 
(e.g. strangeness) production properties close to the phase transition. In 
principle, two distinct event classes (with/without QGP) or the mixed phase 
(coexistence region of hadronic and partonic matter for I s * order phase tran- 
sition) may be reflected in larger event-by-event fluctuations. 

The NA49 experiment used <Jdyn [12] to quantify event-by-event particle 
ratio fluctuations. Odyn measures the difference between widths of parti- 
cle ratio distributions for data and for artificially produced mixed events, 
where only statistical fluctuations are present. Figure [5] (left) shows that 
the dynamical fluctuations in the K/tt ratio are positive. One observes a 
steep rise towards low SPS energies but no significant change from top SPS 
to RHIC energies. The rise towards low energies is not reproduced by the 
UrQMD transport model (there is no significant acceptance dependence). 
The HSD transport model catches the trend but overpredicts the data at 
high SPS energies. Dynamical fluctuations of the K/p ratio (Fig. [5] (right)) 
exhibit two sign changes: a positive plateau at RHIC energies changes to 
a negative plateau at higher SPS energies, followed by a jump to positive 
value at the lowest SPS energy (20,4 GeV). High energy SPS and RHIC data 
are reproduced by the UrQMD model, however the jump at 20 A GeV is not 
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described (this jump between SPS and RHIC is not due to acceptance). The 
values of dynamical event-by-event fluctuations of the p/ir ratio (see con- 
ference slides or [12]) are negative (both at SPS and at RHIC energies) and 
can be reproduced by the UrQMD model (understood in terms of baryon 
resonance decays) in the SPS energy range. 
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Fig. 5. Energy dependence of dynamical K/tt (left) and K/p (right) fluctuations 
for 3.5% most central Pb + Pb (NA49) and for Au + Au interactions (STAR) [12] . 



4. NA49 indications for the Critical Point 



NA61 future data 
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Lattice QCD calculations 
locate the QCD critical point 
at energies accessible at the 
CERN SPS [2]. In "normal" 
liquids (including water) the 
critical point can be rather 
easily detected via the crit- 
ical opalescence phenomenon 
(scattering of light on crit- 
ical long wavelength density 
fluctuations). Over the past 
years several experimental ob- 
servables were proposed to look 
for the CP in heavy ion colli- 
sions. Among them are fluctu- 
ations of mean transverse mo- 
mentum and multiplicity |13j . pion pair (sigma mode) intermittency, ellip- 
tic flow of baryons and mesons, and transverse mass spectra of baryons and 
anti-baryons [H] . One should also note that for strongly interacting matter 
the maximum effect of the CP is expected when the freeze-out happens near 
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Fig. 6. Chemical freeze-out points in NA49 
(squares) and those expected in NA61 (cir- 
cles). See the text for details. 
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the critical point. 

The position of the chemical freeze-out point in the (T — pb) diagram 
can be varied by changing the energy and the size of the colliding system as 
presented in Fig. Therefore we analyzed in NA49 the energy dependence 
of the proposed CP sensitive observables for central Pb+Pb collisions (beam 
energies 2CL4-158A GeV), and their system size dependence (p + p, C + C, 
Si + Si, and Pb + Pb) at the highest SPS energy. Figure [6] shows CP\ 
and CP2 that were considered in NA49 as possible locations of the critical 
point: CP\ with from lattice QCD calculations [2j and T on the empirical 
freeze-out line; CP2 as the chemical freeze-out point of p + p reactions at 
158A GeV (assuming that this freeze-out point may be located on the phase 
transition line). 

4-1. Average pt and multiplicity fluctuations 

At the critical point enlarged fluctuations of multiplicity and mean trans- 
verse momentum are expected [13]. In NA49 we used the & PT fluctuation 
measure [151 ES] an d the scaled variance oj of multiplicity distribution [T7| [T8] 
to study pt and N fluctuations, respectively. For a system of independently 
emitted particles (no inter-particle correlations) <3? PT is equal to zero. For a 
Poisson multiplicity distribution oj equals 1. If A + A reactions are a super- 
position of independent N + N collisions then & PT (A + A) = <£ PT (N + N), 
whereas u(A+A) = u(N+N) + (n) 

^parti where (n) is the mean multiplicity 
of hadrons from ci single N -\- N collisions cmd ojpart 

represents fluctuations 

in N part . The above equations suggest that while <J> PT is independent of 
Npart fluctuations, oj is strongly dependent on them. In the NA49 fixed 
target experiment Np r a ° 3 t can be fixed (spectator energy measured by the 
Forward Calorimeter), whereas cannot be measured. It was shown 

[T9] that fluctuations of iV^f can be suppressed only by selection of very 
central collisions. Therefore multiplicity fluctuations are presented (oj) for 
very central (1%) collisions. 

Figures El El and [TU] present energy (pb) dependence and system size 
(Tchem) dependence of $ PT and oj 0. The chemical freeze-out parameters, 
T c hem{A-, \Jsnn) and ps{A, ^/snn) were taken from fits with the hadron 
gas model [T] to particle yields. The lines correspond to predictions for CP\ 
and CP2 (Fig. [6|) with estimated magnitude of the effects for <E> PT and oj 
at CP\ and CP2 taken from Ref. |13|. [21] assuming correlation lengths £ 
decreasing monotonically with decreasing system size: a) £(Pb+Pb) = 6 
fm and £(p+p) = 2 fm (dashed lines) or b) £(Pb+Pb) = 3 fm and £(p+p) 
= 1 fm (solid lines). The expected magnitudes include NA49 corrections 

3 All $ PT and ui values presented here are obtained in the forward-rapidity region and 
in a limited azimuthal angle acceptance (see corresponding papers for details). 
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due to limited rapidity range (forward-rapidity) and limited azimuthal angle 
acceptance. The width of the enhancement in the (T,hb) plane due to the 
CP is based on Ref. [22J and taken as a(p B ) « 30 MeV and <r(T) ra 10 
MeV. 
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Fig. 7. Energy dependence of $ PT for 7.2% most central Pb + Pb collisions [15] . 
Lines correspond to CP\ predictions (see text); their base-lines are the mean $ PT 
values for 5 energies. 
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Fig. 8. Energy dependence of uj for 1% most central Pb + Pb collisions [18 . Lines 
correspond to CP\ predictions (see text); their base-lines are the mean uj values 
for 5 energies. 

Figures [7] and [8] show no significant energy dependence of pt and multi- 
plicity fluctuations at SPS energies. Thus the results do not provide evidence 
for critical point fluctuations, but a narrower [i B scan would be desirable. 
Figures [9l and [TUl present the system size dependence and show a maximum 
of <£ PT and uj for C + C and Si + Si interactions at the top SPS energy. 
The peak is two times higher for all charged than for negatively charged 
particles as expected for the critical point [T3]. Both figures suggest that 
the NA49 data are consistent with the CP2 predictions. 
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Fig. 9. System size dependence of $ Pr at 158A GeV with p + p, semi-central C + C 
(15.3%) and Si + Si (12.2%), 5% most central Pb + Pb [15]. Lines correspond to 
CP2 predictions (see text) shifted to reproduce the $ PT value for central Pb + Pb 
collisions. 
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Fig. 10. System size dependence of to at 158A GeV for 1% most central p + p |17) . 
C + C and Si + Si [5U] j and Pb + Pb [TB] ■ Lines correspond to CPi predictions 
(see text) shifted to reproduce the w value for central Pb + Pb collisions. 



It was expected that fluctuations due to the CP originate mainly from 
low pt pions [13] . Therefore, in the NA49 analysis of <5 PT the standard 
Pt range (0.005 < pt < 1-5 GeV/c) was divided into two separate pt 
regions: 0.5 < pr < 1.5 GeV/c and 0.005 < pr < 0.5 GeV/c. Indeed, 
the high px region shows fluctuations consistent with zero (see conference 
slides or Ref. [23]) and correlations are observed predominantly at low pr 
(Fig. Hip. However, in the low pt region, data do not show a maximum, but 
a continuous rise towards Pb + Pb collisions. The origin of this behavior is 
currently being analyzed (short range correlations are considered). 
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Fig. 11. The same as Fig. [5] but low px region shown (0.005 < pr < 0.5 GeV/c). 

5. Summary of NA49 results and NA61/SHINE project 

The NA49 experiment obtained numerous interesting results related to 
both the onset of deconfinement and to the critical point. Indications for the 
onset of deconfinement are seen in the energy dependence of the pion yield 
per wounded nucleon N\y (kink), of the {K + ) / (ir + ) ratio and E s (horn), 
and of the mean transverse mass or inverse slope parameters of mx spec- 
tra (step) in central Pb + Pb collisions at lower SPS energies (30 A GeV). 
The results are not reproduced by hadron-string models (RQMD, UrQMD, 
HSD). Extension of the HGM fits the trend of the (K + )/{ir + ) ratio but 
relies on unmeasured hadronic states and educated-guess assumptions on 
the branching ratios. Also the energy dependence of event-by-event hadron 
ratio fluctuations shows interesting effects in the lower SPS energy range 
(increase of adyniK/ir) and sign change of <Jd yn (K/p)), but their relation to 
the onset of deconfinement is not clear. Other (not shown) observables were 
examined for signatures of the onset of deconfinement: a minimum of sound 
velocity at low SPS energies as expected for a 1 st order phase transition [23], 
azimuthal correlations and disappearance of near-side correlations [25] ■ An 
additional energy and system-size scan is important to search for the onset 
of deconfinement in collisions of light nuclei. This is the purpose of the 
NA61/SHINE experiment [26] at the CERN SPS. 

NA49 also searched for indications of the critical point. There are no 
indications of the CP in the energy dependence of multiplicity and mean pj- 
fluctuations in central Pb + Pb collisions. Other (not shown) observables 
were studied: ratios of the anti-baryon/baryon transverse mass spectra and 
elliptic flow V2- Neither shows indications of CP [14] , However, the system 
size dependence of mean pt and multiplicity fluctuations at 158A GeV shows 
a maximum in the complete pt range (consistent with CP2 predictions) and 
an increase from p+p up to Pb + Pb collisions in the low pt region. The low 
Pt region will be carefully analyzed for the effects of short range correlations 
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on <&p T and uj. A detailed energy and system-size scan is necessary to 
establish the existence of the critical point. Therefore, the CP search will 
be continued by the NA61/SHINE experiment. 

The NA61/SHINE experiment (Fig. Q2J) is the successor of NA49. 
The main detector components are inherited from NA49. Several upgrades 
were already completed or are planned in the future: a) in 2007 a forward 
ToF wall was constructed to extend ToF acceptance for particles with p < 
3 GeV/c, b) in 2008 the TPC read-out and Data Aqusition system were 
upgraded to increase the event recording rate by a factor of ~10, c) in 2011 
the NA49 Forward Calorimeter will be replaced by the Projectile Spectator 
Detector (PSD) with a five times better resolution. This excellent resolution 
(about one nucleon in the studied energy range) is crucial especially for the 
analysis of multiplicity fluctuations. 
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In the NA61/SHINE experiment hadron production in p + p, p + A, 
h + A, and A + A reactions at various energies will be analyzed. A broad 
experimental program is planned: search for the critical point, study of the 
properties of the onset of deconfinement, high px physics (energy depen- 
dence of the nuclear modification factor), and analysis of hadron spectra 
for the T2K neutrino experiment and for the Pierre Auger Observatory and 
KASCADE cosmic-ray experiments. Within the NA61/SHINE ion program 
we plan, for the first time in history, to perform a 2D scan with system size 
and energy. The future data onp + p (2009-2010; data taking in progress), 
n B + 12 C (2010-2013), 40 Ar + 40 Ca (2012), and 129 Xe + 139 La (2014) will 
allow to cover a broad range of the phase diagram (see Fig. [6]). With these 
data we will be able to: 



4 SHINE - SPS Heavy Ion and Neutrino Experiment 
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1. Study the properties of the onset of deconfinement. In principle one 
can search for the onset of the 'horn', 'kink', 'step' in collisions of light 
nuclei (the structures observed for Pb + Pb/Au + Au should vanish 
with decreasing system size) 

2. Search for the critical point. An increase of the critical point signal, 
the so-called hill of fluctuations, is expected for systems freezing-out 
near the critical point. Therefore non-monotonic dependence of the 
CP signal on control parameters (energy, centrality, ion size) can help 
to locate the critical point. 

The NA61 program will be complemented by the efforts of other interna- 
tional and national laboratories, BNL RHIC (5 < y / s/v r /v < 39 GeV), JINR 
NICA (3 < < 9 GeV), GSI SIS-100(300) (2.3 < ^sWn < 8.5 GeV) 

and by the heavy ion program at the CERN LHC (y/sNN = 5500 GeV and 
^/s = 14000 GeV for p + p). 
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